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We prove sufficient conditions for Topoiogicai Quantum Order (TQO). The crux of the proof 
hinges on the existence of low-dimensional Gauge-Like Symmetries (GLSs), thus providing a unify- 
ing framework based on a symmetry principle. All known examples of TQO display GLSs. Other 
systems exhibiting such symmetries include Hamiltonians depicting orbital-dependent spin exchange 
and Jahn- Teller effects in transition metal orbital compounds, short-range frustrated Klein spin mod- 
els, and p+ip superconducting arrays. We analyze the physical consequences of GLSs (including 
topological terms and charges) and, most importantly, show the insufficiency of the energy spec- 
trum, (recently defined) entanglement entropy, maximal string correlators, and fractionalization in 
establishing TQO. Duality mappings illustrate that not withstanding the existence of spectral gaps, 
thermal fluctuations can impose restrictions on suggested TQO computing schemes. Our results 
allow us to go beyond standard topological field theories and engineer new systems with TQO. 
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The Landau theory of phase transitions is a landmark 
in physics [l[ . Essential is an order parameter character- 
izing the thermodynamic phases of the system. A new 
paradigm, Topological Quantum Order (TQO), Q ex- 
tends the Landau symmetry-breaking framework. At its 
core, TQO is intuitively associated with insensitivity to 
local perturbations. As such, TQO cannot be described 
by local order parameters. Interest is catalyzed by the 
prospect of fault-tolerant quantum computation Q . 

Several inter-related concepts are typically invoked in 
connection to TQO: symmetry, degeneracy, fractional- 
ization of quantum numbers, maximal string correlations 
(non-local order), among others. The main issue is what 
is needed to have a system with TQO. But a real prob- 
lem is that there is no unambiguous definition of TQO. 
The current article aims to show relations between these 
different concepts, by rigorously defining and establish- 
ing the equivalence between some of them and more lax 
relations amongst others. Most importantly, we (i) prove 
that systems harboring generalized e?-dimensional (with 
d = 0,1,2) Gauge-Like Symmetries (d-GLSs) exhibit 
TQO; (ii) analyze the resulting conservation laws and 
the emergence of topological terms in the action of theo- 
ries in high space dimensions; (iii) affirm that the struc- 
ture of the energy spectrum is irrelevant for the existence 
of TQO (the devil is in the state itself); (iv) establish 
that, fractionalization, string correlators, and entangle- 
ment entropy are insufficient criteria for TQO; (v) report 
on a general algorithm for the construction of string cor- 
relators; (vi) suggest links between TQO and problems 
in graph theory. Our goal is to provide a unifying frame- 
work allowing the creation of new physical models dis- 
playing TQO. A very detailed review and derivation of 
the results presented here is available in 

We focus on quantum lattice systems (and their contin- 
uum extension) having N s = Yl^=i sites, with the 
number of sites along each space direction /i, and D the 



dimensionality of the lattice A. Associated to each lattice 
site (or mode, or bond, etc) i £ Z Ns there is a Hilbert 
space Hi of finite dimension T>. The Hilbert space is the 
tensor product of the local state spaces, H = ^> i Hi, in 
the case of distinguishable subsystems, or a proper sub- 
space in the case of indistinguishable ones. Statements 
about local order, TQO, fractionalization, entanglement, 
etc., are relative to the particular decomposition used to 
describe the physical system. Typically, the most natural 
local language |5( is physically motivated. 

To determine what is needed for TQO, we start by 
defining it. Given a set of N orthonormal ground states 
(GSs) {\g a )}a=i,...,N and a (uniform) gap to excited 
states, TQO exists iff for any bounded operator V with 
compact support (i.e. any quasi-local operator), 



{ga\V\gp) =v 5 al3 + c, 



(1) 



where v is a constant and c is a correction that it is ei- 
ther zero or vanishes exponentially in the thermodynamic 
limit. We will also examine a finite temperature (T > 0) 
extension for the diagonal elements of Eq. (fT]), 

(V) a = tr (p a V) = v + c (independent of a), (2) 

with p a = exp[— H a / (ksT)] a density matrix correspond- 
ing to the Hamiltonian H endowed with an infinitesimal 
symmetry-breaking field favoring order in the state \g a )- 
A system displays finitc-T TQO if it satisfies both Eqs. 
©, and ©. 

A d-GLS of a theory given by H (or action S) is a 
group of symmetry transformations Qd such that the min- 
imal non-empty set of fields 4>i changed by the group 
operations spans a d-dimensional subset C C A. These 
transformations can be expressed as @: t/ik = T7 igCi gik, 
where C\ denotes the subregion 1, and A = (Ji^-i- (The 
extension of this definition to the continuum is straight- 
forward.) Gauge (local) symmetries correspond to d = 0, 
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while in global symmetries the region influenced by the 
symmetry operation is d = D-dimensional. These sym- 
metries may be Abelian or non-Abelian. 

What are the physical consequences of having a sys- 
tem endowed with a symmetry group QdJ Symmetries 
generally imply the existence of conservation laws and 
topological charges with associated continuity equations. 
Wc find that systems with d-GLSs lead to conserva- 
tion laws within d-dimensional regions. To illustrate, 
consider the (continuum) Euclidean Lagrangian density 
of a complex field <f>(x) — (0i(aT), <t>2(x), 03 (x)) (x = 
( Xl ,x 2 ,x s )): £ = |EJ^| 2 + i£Jd T ^| 2 + W(<^), 

with W{^) = |<M 2 ) 2 -i E M ™ 2 (l^l 2 ) and n, v = 

1,2,3. C displays the continuous d = 1 symmetries 
cf)^ — > e i ^^ x "^"^'(f) f j t . The conserved d = 1 Noether cur- 
rents are tensors given by = 5^ - (<9„0* )</> M ], 
which satisfy d = 1 conservation laws [d^j^ + d T j^ T \ = 
(with no summation over repeated indices implicit). 
What is special about d-GLSs is that there is a conser- 
vation law for each line associated with a fixed value of 
all coordinates x u -t^ relating to the topological charge 
Qa»({^m}) = f dx » ^t{x). 

Can we spontaneously break d-GLSs? The absolute 
values of quantities not invariant under Qd are bounded 
from above by the expectation values that they attain in 
a d-dimensional Hamiltonian H (or corresponding action 
S) which is globally invariant under Qd and preserves the 
range of the interactions of the original systems @ . As 
the expectation values of local observables vanish in low-d 
systems, this bound strictly forbids spontaneous symme- 
try breaking (SSB) of non-C/d invariant local quantities 
in systems with interactions of finite range and strength 
whenever d — (Elitzur's theorem), d = 1 for both dis- 
crete and continuous Qd- : and (as a consequence of the 
Mermin-Wagncr-Colcman theorem) whenever d = 2 for 
continuous symmetries 0. Discrete d — 2 symmetries 
may be broken (e.g. the finite-T transition of the D = 2 
Ising model, and the d = 2 Ising GLS of D = 3 or- 
bital compass systems). In the presence of a finite gap 
in a system with continuous d < 2 symmetries, SSB 
is forbidden even at T = [f|. The absence of SSB 
of non-GLSs invariant quantities is due to the presence 
of low-dimensional topological defects such as: domain 
walls/solitons in systems with d = 1 discrete symme- 
tries, vortices in systems with d = 2 U(l) symmetries, 
hedgehogs for d = 2 SU(2) symmetries. Transitions and 
crossovers can only be discerned by quasi-local symme- 
try invariant quantities (e.g. Wilson-like loops) or, by 
probing global topological properties (e.g. percolation in 
lattice gauge theories 0]). Extending the bound of @ 
to T = 0, wc now find that if T = SSB is precluded 
in systems with d-GLSs then it will also be precluded in 
the higher-dimensional system for quantities not invari- 
ant under exact or T = emergent d-GLSs. Exact sym- 
metries refer to [U, H] = 0; in emergent symmetries [H[ 



unitary operators U £ ^emergent are not bona fide sym- 
metries ([[/, H] ^ 0) yet become exact at low energies: 
when applied to any GS, the resultant state must also 
reside in the GS manifold, U\g a ) = E/3 u ai3\9i3)- Thus, 
SSB of T = 0, d = discrete or d < 1 continuous symme- 
tries, in systems of finite interaction strength and range, 
cannot occur. In gapped systems, T = SSB of d < 2 
continuous symmetries is prohibited. 

Degeneracies imply the existence of symmetries which 
effect general unitary transformations within the degen- 
erate manifold and act as the identity operator outside it. 
In this way, GS degeneracies act as exact GLSs. The ex- 
istence of such unitary symmetry operators (generally a 
subset of SU(N)) allows for fractional charge ("iV-ality" 
in the SU (3) terminology of quantum chromodynamics) 
implying that degeneracy allows for fractionalization de- 
fined by the center of the symmetry group. The (m-)rized 
Peierls chains constitute a typical example of a system 
with universal (m— independent) symmetry operators jij, 
where fractional charge quantized in units of e* = e/m 
with e the electronic charge is known to occur The 
bounds above generalize to these symmetries. Different 
Peierls chain GSs break discrete d = 1 symmetries (vio- 
lating Eq. {]]) in this system with fractionalization). The 
Fermi number Nf in the Peierls chain and related Dirac- 
like theories is an integral over spectral functions ; the 
fractional portion of Nf stems from soliton contributions 
and is invariant under local background deformations. 

When the bound of @ is applied anew to correlators 
and spectral functions, it dictates the absence of quasi- 
particle (qp) excitations in many instances [Io| . Here 
we elaborate on this: The bound of 0] mandates that 
the absolute values of non-symmetry invariant correla- 
tors \G\ = | Ofijdlienj 001 with n i c c b and 
c- numbers, arc bounded from above (and from below for 
G > (e.g. that corresponding to (|</>(k, w)| 2 ))) by abso- 
lute values of the same correlators \G\ in a d-dimcnsional 
system defined by Cj. In particular, {a^} can be cho- 
sen to give the Fourier transformed pair-correlation func- 
tions. This leads to stringent bounds on viable qp weights 
and establishes the absence of qp excitations in many 
cases. In high-dimensional systems, retarded correla- 
tors G generally exhibit a resonant (qp) contribution. 
Here, G = G res (k, lo) + G non - re s(k, cj) with G re5 (k, u) = 
u _^_ i0 + ■ bi low-dimensional systems, the qp weight 
Zk — > and the poles of G are often replaced by weaker 
branch cut behavior. If the momentum k lies in a lower 
d-dimcnsional region Cj and if no qp resonant terms 
appear in the corresponding lower-dimensional spectral 
functions in the presence of non-symmetry breaking fields 
then the upper bound 0] on the correlator \G\ (and on 
related qp weights given by lim a ,_ >ek (a; — £k)G(k, cj)) of 
non-symmetry invariant quantities mandates the absence 
of normal qps. Putting all of the pieces together we see 
that if fractionalization occurs in the lower-dimensional 
system then its higher-dimensional realization follows. 
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Our central contention is that in all systems known to 
harbor TQO (and in new examples), d-GLSs are present. 
Old examples include: Quantum Hall systems, Z 2 lattice 
gauge theories, the Toric-code model Q and other sys- 
tems. In all cases of TQO, we may cast known "topo- 
logical" symmetry operators as general low-dimensional 
d < 2 GLSs (e.g. in the Toric code model, there are d = 1 
symmetry operators spanning toric cycles). The presence 
of these symmetries allows for the existence of freely- 
propagating decoupled d-dimensional topological defects 
(or instantons in (d + 1) dimensions of Euclidean space 
time) which eradicate local order. These defects enforce 
TQO. We now state a central result: 

Theorem. When, in a system of finite interaction 
range and strength which satisfies Eq. (Qp, all GSs may 
be linked by discrete d < 1 or by continuous d < 2 GLSs 
U £ Qd, then the system displays finite-T TQO. 

Proof: Let's start by decomposing V = Vq + Vj_. Here, 
Vq is the portion that transforms as a singlet under Qd 
i.e., [U, V ] = 0. To prove the finite-T relation of Eq. ©, 
we write the expectation values over a complete set of 
orthonormal states {|a)}, 



(V }a = lim 



£ a (a|y |a)e 



-f3(E a +<p a a h) 



lim 



lim 



ZMuiv u\i 



j2 e -P(Ea+<t>t, h) 

E b (b\Vo\b)e-^ E "+f b ^ _ 



(Vo)(3. (3) 



Here, we invoked U\a) = \b), and E a — E\, (as [U,H] = 
0) . The term 0° monitors the order parameter content of 
the state \a) vis a vis that preferred by the GS \g a )- In the 



above derivation, 



OJa 



which is evident by 



the application of a simultaneous unitary transformation 
in going from \a) —> u \b) and \g a ) -^ u \gp). V± is not in- 
variant under Qd ([[/, V±\ + 1 0) then by the theorem of 0, 
{y±)a = 0, i.e. for systems with low-dimensional GLSs 
SSB is precluded. Eq. is valid whenever [U, Vq] = 
for any symmetry U. However, [U, Vj_] + 1 implies 
(V±) a = only if U is a low-dimensional GLS. In sys- 
tems in which not all GS pairs can be linked by the use 
of low-dimensional GLSs U £ Qd (U\g a ) = \gp)), finite-T 
SSB may occur. We conclude this proof with two re- 
marks (J]: (i) T = TQO holds whenever all GSs may 
be linked by continuous d < 2 GLSs. (ii) In many systems 
(whether gapped or gapless), T = TQO states may be 
constructed by employing Wigner-Eckart type selection 
rules for GLSs. 

A related consequence of systems with d-GLSs is that 
topological terms which appear in d + 1-dimensional the- 
ories also appear in higher D + 1-dimensional systems 
(D > d). These topological terms appear in actions S 
which bound quantities not invariant under the d-GLSs. 
For instance, in the isotropic D — 2 (or 2+1) dimensional 



general spin tig KK model [i, 11 [, of exchange constant 
J > 0, the corresponding continuum Euclidean action is 

of the 1+1 form S = ^- J dxdr ^-{d T m) 2 —v s {d x mf 



(d^m-d^m), with v € {x, r} and 



the rank two Levi Civita symbol. Here, as in the non- 
linear-CT model of a spin-5 chain, m a normalized slowly 
varying staggered field, g = 2/5, v s = 2JS, 9 = 2ttS, and 
Str a "transverse-field" action term which does not act on 
the spin degrees of freedom along a given chain Q is 
the Pontryagin index corresponding to the mapping be- 
tween the 1+1 dimensional space-time (x,t) plane and 
the two-sphere on which rh resides. This 1+1 dimen- 
sional topological term appears in the 2+1 dimensional 
KK system even for arbitrary large positive coupling J. 
This, in turn, places bounds on the spin correlations and 
implies, for instance, that in D = 2 integer-spin t2 g KK 
systems, a finite correlation length exists. 

What characterizes TQO? We will see that neither the 
spectrum of H, nor the entanglement entropy [l^ |. nor 
string correlations are sufficient criteria. That the spec- 
trum, on its own, is insufficient is established by counter- 
examples, e.g. that of the D = 3 Z 2 lattice gauge theory 
which is dual to a D = 3 Ising model [l]|. Albeit shar- 
ing the same energy spectrum, the gauge theory displays 
TQO while the Ising model harbors a local order. A spec- 
tral equivalence exists between other TQO inequivalent 
system pairs. For example, Kitaev's Q and Wen's models 



14j are equivalent and have a spectrum which is identi- 



cal to that of an Ising spin chain with nearest neighbor 
interactions [3] • The mappings between these systems to 
an Ising chain demonstrate that despite the spectral gap 
in these systems, the toric operator expectation values 
may vanish once thermal fluctuations arc present. These 
mappings also illustrate that the quantum states them- 
selves in a particular (operator language) representation 
encode TQO and that the duality mappings, being non- 
local in the original representation, disentangle the order. 
Thus, we cannot, as coined by Kac, "hear the shape of a 
drum". The information is in the eigenvectors. 



It was recently suggested |12| that a non-local borne 
deviation 7 > from an asymptotic area law scaling for 
the entanglement entropy constitutes a sharp measure of 
TQO. We find, for example in Klein spin models which 
do not host TQO (i.e. do not satisfy Eq. (JTJ) for the 
set of all GSs), an arbitrary - contour shape dependent 

entanglement entropy which deviates from an area law 
[3| • This suggests that the criterion of can be violated 
and may require additional improvements, such as the 
specification of an exact limit and/or average for large 
contours. 

The insufficiencies of the spectra and entropy in deter- 
mining if TQO is present have counterparts in the topol- 
ogy of graphs and in the Graph Equivalence Problem 
(GEP) in particular. The adjacency matrix of a graph 
has elements Cy = 1 if vertices i and j are linked by an 
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edge and Gy = otherwise. Vertex relabeling i — ► p(i) 
leaves a graph invariant but changes the adjacency ma- 
trix G according to C — ► C — P^CP with P an orthogo- 
nal matrix which represents the permutation: P = <5jjp(i)- 
The GEP is the following 0: "Given C and G', can we 
decide if both correspond to the same topological graph?" 
The spectra of G and C are insufficient criteria. Entropic 
measures 



151 are useful but also do not suffice. 



In many systems (with or without TQO), there are 
non-local "string" correlators which display enhanced (or 
maximal) correlations vis a vis standard 2-point correla- 
tion functions. We now outline an algorithm for the con- 
struction of such non-local correlators. [In systems with 
uniform global order already present in their GS, the al- 
gorithm leads to the usual 2-point correlators.] We seek 
a unitary transformation U s which rotates the GSs into a 
new set of states which have greater correlations as mea- 
sured by a set of local operators {VI}. These new states 
may have an appropriately defined polarization (eigenval- 
ues {i>i}) of either (i) more slowly decaying (algebraic or 
other) correlations, (ii) a uniform sign {partial polariza- 
tion) or (iii) maximal expectation values Vi = v max for all 
i (maximal polarization). Cases (i) or (ii) may lead to a 
lower dimensional gauge like structure for the enhanced 
correlator. In systems with entangled GSs (such as in 
many (yet not all) GSs hosting T = TQO), U s cannot 
be a uniform product of locals; U s ^ JlieA ^i- To provide 
a concise known example where these concepts become 
clear, we focus below on case (ii) within the well-studied 
AKLT Hamiltonian [l6[. Here, there is a non-trivial uni- 
tary operator U s = I^eaq^iTrSfSft, {[U„H] ? 0) 
which maps the GSs {|<7ct)} into linear superpositions of 
states in each of which the local staggered magnetiza- 
tion VI = (— iyS? = M\ is uniformly non- negative (or 
non-positive) at every site i. As all transformed states 
\Pa) = U s \g a ) are superpositions of states with uniform 
sign Vi (allowing for two non-negative or non-positive v\ 
values at every site out of the three 5=1 states) , Gy = 



(PalViVjK) > (<7 a |V]Vj|<7„) ee G K 
ten in full, Gy becomes a non-local 
Gy = (-l)^(g a \S? ni< k<j expK^?| 9Q 
for arbitrarily large separation |i — j|. In the 
basis, there are 2 N ° states \<j>) with non-negative (non- 
positive) Vi for all i (they form a Hilbert subspace), and 



When U s is writ- 

"string" correlator 

> = (2/3) 2 
®S? eigen- 



all of these have the same value of G, 



\ViVi\ 



that 



is 4/9. Such an exponential number of states appears in 
systems with a local Z 2 gauge structure. Here, the string 
correlators exhibit a symmetry under local (d = 0) gauge 
transformations while H itself is not invariant. These 
d = Z2 transformations correspond (in the ®S( ba- 
sis) to the creation (annihilation) of an S? = state 
at any site j followed by a unit displacement of 5^ at 
all sites k > j. The operators linking the GSs form 
a discrete (Z 2 x Z 2 ) global symmetry group [lj| and, 
as such, T = SSB may occur. The general polariza- 



tion operators are intimately tied to selection rules which 
hold for all of the GSs; these selection rules (and a par- 
ticular low-energy projection of truncated Hamiltonians 
for gapped systems) enable the construction of general 
string and higher dimensional "brane" correlators. In the 
AKLT problem, the local expectation value (g a \Sf =1 \g a ) 
depends on \g a ), violating Eq. JTJ) and suggesting that 
the AKLT chain is not topologically ordered. We may 
also generate maximal (Gy = 1) string correlators; here, 
the number of states with maximal polarization is finite 
and a non-local gauge-like structure emerges. Generally, 
in any system with known (or engineered) GSs, we may 
explicitly construct polarizing transformations U s . 

To conclude, we proved a theorem establishing how 
d < 2 dimensional GLSs in gapped systems can man- 
date TQO via freely propagating low-dimensional topo- 
logical defects (e.g. d = 1 solitons). This result extends 
the set of known systems which exhibit TQO to include 
new orbital, spin, and Josephson junction array prob- 
lems. All known TQO systems display GLSs We ex- 
amined physical consequences of GLSs (such as conserva- 
tion laws and topological terms in high-Z? theories), and 
low-dimensional gauge-like structures in theories with en- 
tangled GSs in non-maximal string correlators (irrespec- 
tive of TQO). We established insufficiencies of the Hamil- 
tonian spectra, fractionalization, entanglement entropy, 
and maximal string correlations in determining TQO. 
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